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Abstract We present the results of an extensive study of the detectability of Earth- 
^ . sized planets and super-Earths in the habitable zones of cool and low-mass stars using 

transit timing variation method. We have considered a system consisting of a star, 
■^r . a transiting giant planet, and a terrestrial-class perturber, and calculated TTVs for 

different values of the parameters of the system. To identify ranges of the parameters 
for which these variations would be detectable by Kepler, we considered the analysis 
P I ■ presented by Ford et al. (2011, I arXiv: 1102.0544 ) and assumed that a peak-to-peak 

^n \ variation of 20 seconds would be within the range of the photometric sensitivity of 

this telescope. We carried out simulations for resonant and non-resonant orbits, and 
identified ranges of the semimajor axes and eccentricities of the transiting and per- 
turbing bodies for which an Earth-sized planet or a super-Earth in the habitable zone 
O ' of a low-mass star would produce such TTVs. Results of our simulations indicate that 

in general, outer perturbers near first- and second-order resonances show a higher 
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^ ■ prospect for detection. Inner perturbers are potentially detectable only when near 1:2 

and 1:3 mean- motion resonances. For a typical M star with a Jupiter-mass transiting 
planet, for instance, an Earth-mass perturber in the habitable zone can produce de- 
tectable TTVs when the orbit of the transiting planet is between 15 and 80 days. We 
present the details of our simulations and discuss the implication of the results for the 
detection of terrestrial planets around different low-mass stars. 
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1 Introduction 



Cool and low-mass stars (such as M dwarfs) present the most promising candidates 
- , , , for searching for Earth-like planets and super-Earths. These stars show large reflex 
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accelerations in response to an orbiting body, and their low surface temperatures, 
which place their liquid water habitable zones (Kasting et al 1993) at close distances 
of 0.1-0.2 AU (corresponding to orbital periods of 20-50 days), cause their light-curves 
to show large decrease when they are transited by a close-in planet. In the past few 
years, these advantages have resulted in the detection of several super-Earths around 
low-mass stars including a 6.0-7.5 Earth-mass planet around the star GJ 876 (Rivera 
et al. 2005, 2010; Correia 2010), 4-6 planets with masses ranging from 1.7 to 7.0 Earth- 
masses around star GL 581 (Bonfils et al 2005; Udry et al 2007; Mayor et al 2009; Vogt 
et al 2010), and two transiting planets around M stars GJ 436 (Butler et al 2004; Gillon 
2007) and GJ 1214 (Charbonneau 2009). 

The existence of systems such as GJ 876 and GL 581, where an M star is host 
to multiple planets in short-period orbits, suggests that other single-planet low-mass 
stars may also host additional objects. In systems where one (or more) of these bodies 
transits, the interaction between this object and other planets of the system results in 
the exchange of energy and angular momentum, and creates small, short-term oscilla- 
tions in their semimajor axes and orbital eccentricities. On the transiting planet, these 
oscillations manifest themselves as variations in the time and duration of the planet's 
transit. In this paper, we discuss the possibility of the detection of these transit timing 
variations when they are generated by an Earth-like planet or a super-Earth in the 
habitable zone of a low-mass star. 

Measurements of transit timing variations can be used in detecting small, non- 
transiting planets (Miralda-Escude 2002, Schneider 2003, Holman and Murray 2005, 
Agol et al 2005, Heyl and Gladman 2007). Despite the complication in extracting the 
mass and orbital elements of the perturbing body from the values of TTVs (Nesvorny 
and Morbidelli 2008, Nesvorny 2009, Nesvorny and Beauge 2010, Meschiari and Laugh- 
lin 2010, Veras et al 2011), this method has also been recognized as a potential mech- 
anism for detecting Trojan planets (Ford and Gaudi 2006, Ford and Holman 2007, 
Haghighipour et al, in preparation), as well as the moons of transiting gas-giants 
(Sartoretti and Schneider 1999, Simon et al 2007, Kipping 2009a&b). TTVs are also 
sensitive to the precession (Miralda-Escude 2002) and inclination of the orbit of the 
transiting object (Payne et al 2010), and in close, eclipsing binary systems, they can 
be used to detect circumbinary planets (Schneider and Chevreton 1990, Schneider and 
Doyle 1995, Doyle et al 1998, Doyle and Deeg 2003, Deeg et al 2008, Sybilski et al 
2010, Schwarz et al 2011). 

In this paper, we focus our study on the detection of TTVs using Kepler space 
telescope. Kepler monitors more than 150,000 stars with spectral types of A to M, 
and with apparent magnitudes brighter than 14 (Oshagh et al 2010; Batalha et al 
2010, Coughlin et al 2011). A survey of the discovery- space of this telescope indicates 
that low-mass stars such as M dwarfs, constitute a significant number of its targets. 
The large number of the targets of this telescope, combined with its unprecedented 
photometric precision, and the long duration of its operation, makes Kepler ideal for 
searching for transit timing variations. A recent success, for instance, has come in the 
form of detecting variations in the orbital periods of transiting planets around stars 
Kepler-9 (Holman et al 2010) and Kepler- 11 (Lissauer et al 2011). Further analysis of 
the currently available data from this telescope also points to the existence of many 
additional candidates that have potentially detectable variations in the times of their 
transits (Steffen et al 2010, Ford et al 2011). 

As shown by Agol et al (2005), Steffen and Agol (2005), and Agol and Steffen 
(2007), the amplitude of TTV is greatly enhanced when the transiting and perturbing 



planets are near a mean- motion resonance (MMR). In this case, transit times vary as 
a function of the ratio of the masses of the two planets which enhances the sensitivity 
of the TTV method by one to two orders of magnitude. For instance, as shown by 
Agol et al (2005), around a solar-mass star, an Earth-mass planet can create TTVs 
with amplitudes of up to 180 s when near an exterior 2:1 MMR with a transiting 
Jupiter-size body in a 3-day orbit. In the recently detected Kepler-9 system where two 
Saturn-sized planets are near a 2:1 resonance, the rate of the transit timing variation on 
the inner planet was even higher and reached to approximately 50 minutes. This is two 
orders of magnitude higher than the '-^lO s transit timing precisions that was achieved 
by a 1.2 m ground-based telescope in the measurements of the transit timing of HD 
209458 (Brown et al 2001). The fact that a perturber near a mean-motion resonance 
can produce large TTVs implies that such resonant systems will have a higher chance 
of being detected by the TTV method. In this paper, we will pay especial attention to 
such resonant orbits. 

The outline of this paper is as follows. In section 2, we present the details of our 
methodology for calculating TTVs. Since we are interested in the detectability of an 
Earth-like planet or a super-Earth in the habitable zone of a star, we will analyze the 
dependence of the magnitudes of TTVs on different parameters of a system and present 
the results in section 3. Section 4 has to do with the implications of this analysis for 
the detection of planets in the habitable zone, and section 5 concludes this study by 
presenting a summary and reviewing the results. 



2 Calculation of TTV 

The system of our interest consists of a star, a transiting planet, and a perturbing 
body. Since the focus of our study is on low-mass stars, we choose the star to have 
a mass between O.25M0 and O.TM©. The transiting planet is taken to be initially on 
a circular orbit with a period (Ptr) between 3 and 10 days. The mass of this planet 
(mtr) is varied between 0.1 and 3 Jupiter-masses (Mj). The perturber is considered to 
be smaller with a mass (?7ipr) ranging from 0.3 to 10 Earth-masses (Mq). The initial 
orbital eccentricity of this planet is varied from to 0.3 in increments of 0.05, and its 
orbital period (Ppr) is chosen such that 0.3 < Ppr/Ptr < 10. We consider the system to 
be coplanar, and unless stated otherwise, we set the initial values of the mean-anomaly, 
argument of pericenter, and the longitude of the ascending node of the transiting planet 
equal to zero. 

We integrated our systems using the Bulirsch-Stoer algorithm of the N-body in- 
tegration package Mercury (Chambers 1999). This algorithm was modified to include 
general relativistic effects (Saha and Tremaine 1992). We carried out our integrations 
for three time spans of 1, 3, and 30 years. The 3-year time span was chosen to correspond 
to the duration of Kepler's primary mission. The 1-year and 30-year integrations were 
meant to show the short-term and long-term effects. The 30-year time span was chosen 
to ensure that long-term variations in the amplitude of TTVs will also be modeled. 

To calculate TTVs, we obtained the difference between the time of the mid-transit 
in the unperturbed system (star and transiting planet), ti, and its corresponding value 
obtained from the numerical integrations of the perturbed system (star, transiting and 
perturbing planet), t2. We assumed that at t = 0, the centers of the star and the 
transiting planet were on the x-axis and calculated the time of transit by interpolating 
between the times before and after the center of transiting planet crossed the center 



Table 1 Initial orbital elements of the perturber in figurel. 



Orbital Elements Top Middle Bottom 



Semimajor Axis (AU) 


0.089 


0.089 


0.042 


Eccentricity 











Inclination 











Argument of Periastron 


352.8° 


352.8° 





Mean- Anomaly 


7.192° 


7.192° 





Longitude of Periastron 


352.8° 


352.8° 





Longitude of Ascending Node 












of the star. Figure 1 shows a sample of our results. The system in this figure consists 
of a O.367M0 star, a Jovian- type transiting planet with a mass of 1.05Mj on a 10-day 
orbit, and an Earth-mass perturber. Table 1 shows the orbital elements of this object. 

The top two panels in figure 1 show the TTVs for a non-resonant case where the 
perturber is in an exterior orbit with a period of Ppr — 16.11 days. As expected, 
the amplitude of TTVs in this system is small. The results shown here are consistent 
with the analytical prediction of TTVs as given by Agol et al (2005). These authors 
indicated that in a non-resonant system, when the orbital period of the perturber 
is half-way between two resonances, the timing deviation is approximately given by 
0.7Ptr(^^pr/M)[atr/(ttpr — Citr)] where M is the mass of the central star. For the 
top two panels of figure 1, the orbital period of the perturbing planet (16.11 days) is 
between the 4:3 and 2:1 exterior resonances. In this case, the formula above predicts 
a timing deviation of approximately 37 s which is consistent with the RMS in these 
two panels. The bottom panel of figure 1 corresponds to a case of near interior 1:2 
mean-motion resonance. As expected, the amplitude of TTVs in this case vary in a 
large range with a maximum of ~800 s. 

As shown by Agol et al (2005), high-order resonances result in smaller magnifica- 
tions in the amplitude of TTV and as such they are more suitable for demonstrating 
the contrast between the short- and long-term variations in transit timing. To show 
this effect, we placed the perturber of the system of figure 1 near an exterior 3:1 
resonance and integrated the system for a longer time. Table 2 shows the orbital el- 
ements of the perturber at the beginning of the integrations. The results are shown 
in figure 2. The top panel in this figure shows TTVs for 250 days. As shown here, 
although the planets are near a resonance, the amplitudes of TTVs are much smaller 
than those of the 1:2 MMR (bottom panel of figure 1). This is an expected result as 
the 3:1 MMR is one order higher than the 1:2 MMR. As shown by Agol et al (2005), 
the analytically predicted value of the amplitude of TTVs in this case is given by 
(Ptr/2 / 7r)(mpr/M)(apr/atr) '^ 1.67 s which agrees with the amplitude of TTVs as 
shown in this graph. 

The continuation of integrations to 1000 days (middle panel) reveals the onset 
of long-term secular variations. As the integrations continue to 10 days, the secular 
long-term variations are fully developed and reach an amplitude as high as 253 s. The 
increase in the amplitude of TTVs in this case can be attributed to the increase in the 
eccentricity of the perturber as a result of its long-term interaction with the massive 
transiting planet. Figure 3 show the eccentricities of these objects for the duration of 
integration. 

Figure 4 shows the maximum values of the amplitudes of TTVs for the system 
of figure 1 (a O.367M0 star, a Jovian- type transiting planet with a mass of 1.05Mj 
on a 10-day orbit, and an Earth-mass perturber) when the perturber is in different 



Table 2 Initial orbital elements of the perturber in figure 2. 



Orbital Elements 


All Panels 


Semimajor Axis (AU) 

Eccentricity 

Inclination 


0.137 




Argument of Periastron 


354.81° 


Mean- Anomaly 
Longitude of Periastron 
Longitude of Ascending Node 


5.194° 

354.81° 





resonant and non-resonant orbit qj. As shown here, low-order resonances show high 
variations in transit timing that are well within the range of the photometric precision 
of Kepler space telescope. These results are also supported by the analytical formulae 
given by Agol et al (2005). These authors have shown that when the two planets 
are in a first order j : j -\- 1 MMR, the maximum value of the TTVs is approximately 
given by (P/4.5j)(mpert/''^pert + ^tr)- For the system of figure 1 when the two planets 
are in an exterior 2:1 resonance, the maximum TTV as predicted by this formula is 
approximately 600 s which agrees with the value shown in figure 4. More results can 
be found in the appendix where similar simulations have been carried out for resonant 
and non-resonant orbits in a system with a Sun-like star, and for a perturber with a 
mass of 1 M0. 



3 Dependence on the parameters of the system 

The magnitude of TTV changes with the mass and orbital elements of the planets, 
and the type of the central star. To identify ranges of these parameters that produce 
detectable TTVs, we carried out simulations for different values of the mass and orbital 
eccentricity of the perturber, period and mass of the transiting planet, and the mass 
of the host star. Calculations were carried out for 3 years and for both near resonance 
and non-resonant orbits. Except when studying the effect of the orbital eccentricity 
of the perturber, calculations began with both planets in circular orbits. We present 
the results of our simulations in the following. Except for the case where the mass of 
the central star is varied, we use the system of figure 1 (?Tistar = O.367M0 , mtr = 
IMj , TTipr = 1^© , Ptr = 10 days) as a sample system for our analysis. 



3.1 Mass of the perturbing planet (mpr) 

Since the interaction between the two planets is gravitational, the mass of the per- 
turber plays an important role in varying the time of transits. We used the system 
mentioned above and calculated TTVs for different values of this parameter. Since we 
are interested in Earth-like planets and super-Earths, we considered different values for 



^ We did not show the case of 1:1 MMR in figure 4 because our integrations indicated that 
in a large area around this resonance, the orbit of the perturbing planet was unstable. We 
carried out a more detailed study of this resonance and were able to identify ranges of orbital 
elements for which a perturber near a 1:1 MMR would have a stable orbit. We will present 
the analysis of transit timing variations for this system, and the detectability of their TTVs, 
in an up-coming article (Haghighipour et al. in preparation). 



rripr between O.IM0 to IOMq. Figure 5 shows the results for two cases of near exterior 
resonances (2:1 and 5:2 MMRs), and for a non-resonant system. As shown here, near a 
resonance, variations in the time of transit scale with the mass of the perturber. This 
is an expected result that has also been reported by Agol et al (2005). Our simulations 
suggest that such a scaling seems to exist for non-resonant orbits as well, which is a 
natural consequence of the conservation of energy and momentum (see equation 30 of 
Agol et al 2005). 



3.2 Eccentricity of the perturbing planet (epr) 

As explained in the beginning of this section, we considered the initial orbits of the 
perturbing and transiting planets in all our simulations to be circular. During the 
course of an integration, the mutual interactions of the two planets cause their orbital 
eccentricities to increase (e.g., figure 4). In cases where the planet-planet interaction 
does not make the system unstable, the increase in the orbital eccentricity of the 
perturbing body results in a decrease in its periastron distance which in turn causes 
rapid and fast changes in the time of transit. To examine the enhancing effect of epr 
on TTVs and its role in their possible detection, we carried out simulations where the 
initial orbital eccentricity of the perturbing planet was chosen to be between and 
0.3. We considered a system with a Jupiter-mass transiting planet in a 10-day orbit, 
an Earth-mass perturber, and a 0.367 Mq central star. We chose the system to be 
co-planar, and set all orbital angles for both planets equal to zero. Results indicated 
that in general, an increase in perturber's orbital eccentricity increases the amplitude 
of TTV. However the rate of this increase depends on the resonant state of the system 
(figures 6 and 7) , and the time of integration (figure 8) . 

For large values of the perturber's eccentricity, high-order MMRs play an important 
role in the variations of transit timing (Pan and Sari 2004; Agol et al. 2005). This can 
be seen in figure 6 where the amplitudes of TTVs are shown for different values of the 
ratios of the orbital periods of the two planets. As shown here, among the resonant 
cases, the rate of the increase in TTVs is larger for a near 5:2 MMR. Figure 6 also shows 
that when the two planets are near a MMR, a range of the perturber's eccentricity exists 
for which the increase in the value of TTV is negligibly small. This can be seen in figure 
7 where the magnitudes of TTV and the period-ratio are shown for a near exterior 2:1 
MMR. As shown here, the amplitude of TTV maintains an almost constant maximum 
value around 600 s for epr < 0.08 (for the system considered in these simulations, 
integration for larger values of eccentricity resulted in instability) . Figures 6 and 7 are 
also consistent with figure 5 of Agol et al (2005) where these authors show the dispersion 
in timing variations for different values of the outer planet's orbital eccentricity and 
for different resonances. Calculations by Agol et al (2005) indicate that when the 
perturber is in a low-eccentricity resonant orbit, the variations of its eccentricity that 
are caused because of its interactions with the massive transiting planet, are of the 
order {rripert/'mtr) ' . For the system of figure 7, this quantity is approximately 0.14. 
Therefore, for initial eccentricities up to half of this value (i.e., ^ 0.07, as in figure 7), 
the TTVs are primarily dominated by the period variations of the perturber and not 
strongly affected by its initial eccentricity. 

To explore how the amplitude of TTVs vary as the time of integration increases in 
a system where the perturber is in a resonant and eccentric orbit, we calculated TTVs 
for different values of the orbital eccentricity of an Earth-mass planet in a 1:2 MMR 



with a Jupiter-mass transiting body. The mass of the central star is 0.367 Mq. Figure 
8 shows the results for different integration times (in terms of the number of transits). 
As shown here, the strongest increase in TTVs is reached when the time of integration 
varies between 20 and 30 transits. This can be explained noting that for a system in 
a j : j + 1 MMR, the libration period scales as 0.5 j~ / (^pert/''^tr)~ ^ (Agol et al 
2005). For the system of figure 8 where j — 1, this results in approximately 24 transits. 
For longer times of integration, the magnitude of TTV increases almost linearly with 
the perturber's eccentricity. 



3.3 Period of the transiting planet (Ptr) 

We calculated TTVs for different values of the period of the transiting planet. Figure 
9 shows a sample of results for a non-resonant (top) and a resonant (bottom) system. 
Simulations were carried out for a sub-Earth perturber and several super-Earths. As 
shown here, the magnitude of TTV increases with increasing the period of the transiting 
planet. This result agrees with the analytical analysis by Holman and Murray (2005) 
and Agol et al (2005) where these authors have shown that near a MMR, the magnitude 
of TTVs scale as Ptr- We will explain the implications of these results for the detection 
of habitable planets in the next section. 



3.4 Mass of the transiting planet (mtr) 

Simulations were also carried out for transiting planets with masses ranging from O.lMj 
to 3Mj. Results point to different mtr-dependence based on the resonant state of the 
system and the order of the resonance. Integrations indicated that for the system 
considered here (i.e. the system of figure 1), except for the interior 1:2 and exterior 
2:1 MMRs, other first-order resonances were unstable. Figure 10 shows the results for 
when the system is near these two resonances. As shown here, the amplitude of TTV 
is larger for smaller values of the mass of the transiting planet. As the mass of this 
object increases, the corresponding value of its TTV signal becomes smaller. This is 
a result that agrees with the theoretical analysis of Agol et al (2005). As shown by 
these authors, the maximum variation of transit timing in a resonant case scales as 
mpr/(mpr + mtr). 

The top panel of figure 11 shows the values of TTV when the system is near two 
high-order resonances of 3:1 and 5:2, and for some of their neighboring non-resonant 
orbits. As shown here, the values of TTVs are slightly increased for large transiting 
planets. This can be attributed to the fact that as the mass of the transiting planet 
increase, its interaction with the perturbing body excites the eccentricity of this object 
which in turn causes higher order resonances to become important. The increase in 
the eccentricity of the perturbing planet can be seen in the bottom panel of figure 11 
where the median of the eccentricity of this object has been plotted against the mass 
of the transiting planet for a sample non-resonant system and also when the system is 
in a 5:2 MMR. A shown here, the increase in the median eccentricity is larger in the 
non-resonant case indicating that, as shown in the top panel, TTVs would be larger 
outside resonances. The decreasing trend of TTV in the case of 5:2 MMR may also be 
due to similar effects as in the case of 2:1 MMR (figure 10). 



3.5 Mass of the central star 

Although the focus of our study has been on the possibility of the detection of TTVs 
around low-mass stars, we carried out simulations for different values of the mass of 
the central star as well. Figure 12 shows the results for a perturber with different 
masses. The graph on the top is for an exterior 2:1 MMR. The graph on the bottom 
corresponds to a non-resonant system where the perturber is in a 28.3-day orbit. As 
shown here, in the non-resonant system, the magnitude of TTV seems to be inversely 
proportional to the mass of the central star. This result has also been reported by 
Agol et al (2005). In the resonant system, on the other hand, the mass of the central 
star does not seem to play an important role in changing the magnitude of TTV. This 
characteristic of resonant systems enables us to examine the effect of the mass of the 
star on the possibility of the detection of the perturber in its habitable zone. We explain 
this in more detail in the next section. 



4 Implications for the detection of a perturber in the habitable zone 

Using the results of the simulations of section 3, we assessed the detect ability of TTVs 
that were produced by a perturber in the habitable zone. We determined the locations 
of the boundaries of the habitable zone (r) using the equation 

where L is the luminosity of the central star. The quantity Vs in this equation, cor- 
responds to the (inner, outer) boundaries of the Sun's habitable zone. In general, this 
quantity is model dependent and varies with different atmospheric circulation models. 
For instance, as suggested by Kasting et al (1993), rs = (0.95 , 1.37) (AU), whereas 
models by Menou and Tabachnik (2003) and Jones et al (2005, 2006) suggest an inner 
boundary as close as 0.7 AU and an outer boundary beyond 2 AU. In this study, we 
considered some intermediate values, and assumed that rs = (0.78 , 1.8) (AU). 

To calculate the luminosity of the star, we used two different mass-luminosity rela- 
tions. First we considered L ^ M ' , where both the luminosity and mass are in solar 
units. We also used the mass- luminosity relation suggested by Duric (2004): 

- L{Lq) - 0.23 M^-^ (Mo) , M < 0.4 M© 

- L{Lq) - M^ (Mo) , M > O.4M0 . 

Calculations indicate that the locations of the boundaries of HZ, as obtained from 
these two relations, are very close to one another. We, therefore, used the boundaries 
obtained from the L ^ M ' for the purpose of our analysis. 

Figures 13 and 14 show the HZ for three stars with masses 0.252, 0.367, and 
O.486M0. Since we are interested in identifying detectable TTVs in systems where 
the transiting planet is at different orbital distances, and also because resonant orbits 
create larger TTVs, we chose the horizontal axis to represent the period of the tran- 
siting body and identified different mean-motion resonances. The shaded area in each 
graph corresponds to systems where the variations of TTVs due to a perturber in an 



initially circular orbit are larger than 20 s. This value is consistent with the photometric 
precision of Kepler as reported by Ford et al (2011). 

An interesting result shown in figures 13 and 14 is the large contrast between the 
sizes of the regions where inner and outer perturbers produce detectable TTVs. This 
can be seen for an Earth-mass perturber (figure 13 and the top panel of figure 14) as 
well as a IOM0 super-Earth (bottom panel of figure 14). As shown in these figures, 
outer perturbers seem to have higher potential for detection. An inner perturber also 
produces detectable TTVs. However, as expected, the orbit of the planet in this case is 
primarily limited to the regions near first- and second-order mean-motion resonances 
(1:2 and 1:3 MMRs). The perturber can also produce detectable TTVs when in the 
habitable zone of the central star. For instance, as shown by the top panel of figure 
13, a IM0 perturber can produces detectable TTVs when in the HZ of a O.252M0 
star and near a 1:3 MMR with a IMj transiting planet with a period larger than 40 
days. As the mass of the star increases, the detect ability in the HZ is shifted to farther 
orbits. 

As expected, an outer perturber produces detectable TTVs in a larger range of 
distances. The region of detect ability in this case contains many first- and high-order 
MMRs. We have shown some of these resonances in figures 13 and 14. As shown here, 
the lower boundary is limited to the exterior 3:2 MMR whereas the upper boundary is 
slightly above the 2:1 MMR and in the case of larger perturbers it may even reach to 
the 3:1 MMR (bottom panel of figure 14). As shown in these figures, a small planet in 
the HZ of a typical M star can produce detectable TTVs on a transiting Jupiter-mass 
body in an orbit with a period larger than 15 (10) days corresponding to an Earth-mass 
(IOM0 super-Earth) perturber. 



5 Summary and concluding remarks 

We carried out an extensive study of the possibility of the detection of small planets 
in the HZ of low-mass stars with Kepler space telescope and transit timing variation 
method. With its high precision photometry, Kepler has the capability of detecting 
small TTVs such as those created by terrestrial-class perturbers. We assumed that 
TTVs with amplitudes as small as 20 s can be detected by Kepler, and identified 
regions of the parameter space where the variations in the time of the transit of a giant 
planet in a short-period orbit around a low-mass star would be larger than this value. 
In general, the lowest detectable amplitude of TTVs depends upon several quanti- 
ties including the magnitude of the central star, the depth of the transit, the durations 
of ingress and egress, and the number of transits observed. There will also be sys- 
tematics that will affect the sensitivity of the telescopes at different times. The 20 s 
detectability limit considered in this study corresponds to the best transit timing pre- 
cision obtained by Ford et al (2011) in their statistical analysis of the transit timing of 
Kepler^s planetary candidates. Analyzing the data obtained from 1235 potential tran- 
siting planets (Borucki et al 2011) during the first four months of the operation of this 
telescope. Ford et al obtained a transit timing accuracy ranging from approximately 
20 s to 100 minutes (figure 1 of Ford et al 2011). The 20 s precision is also consistent 
with the analysis of the light curves of Kepler-4 through Kepler-8 as given by Kipping 
and Bakos (2011). As shown by these authors, the timing precision for the transiting 
planets in these systems ranges between 20 and 250 seconds. It is worth noting that 
systems with TTVs larger than 15-20 s may also be detected by CoRoT. This telescope 
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has been monitoring more than 12000 dwarf stars with apparent magnitudes ranging 
from 11 to 16. 

Although many non-resonant systems produced detectable TTVs, the results of 
our study indicate that as expected, perturbers near mean-motion resonances with 
the transiting planet show a higher prospect of detection. This can also be seen in 
the system of Kepler-9 where the two TTV-detected planets of this system are near a 
2:1 MMR (Holman et al 2010). Our study shows that interior perturbers with masses 
ranging from 1 to 10 Mq can produce detectable TTVs when near 1:2 and 1:3 MMRs. 
However, in the majority of our systems, these interior perturbers were not in the HZ. 
When in an outer orbit, these planets can be in the HZ and produce detectable TTVs 
when in first- and second-order resonances. For instance, for a typical M star with a 
mass of approximately 1/3 of the Sun where the HZ extends between 0.1 and 0.2 AU, 
perturbers in 3:2, 2:1, 5:3, and 3:1 MMRs produced large-magnitude TTVs when in 
the system's HZ. For the systems studied here, a perturber with a mass ranging from 
1 to 10 Earth-masses can be detected in the HZ when the transiting planet is in 15-80 
day orbits. 

In this study we considered the planetary system to be planar and consists of only 
three bodies. These assumptions may limit the applicability of our results. For instance, 
an inclined perturber will cause the inclination of the orbit of the transiting planet to 
vary which in turn affects the variations in the times of its transits (Payne et al 2010). 
Also, as suggested by the simulations of planet formation around low-mass stars, these 
stars may be hosts to more than two planets (e.g., GJ 876 hosts four planets, and GL 
581 has 4-6 super-Earths). Depending on their masses and orbital elements, additional 
bodies may play a role in the magnitude of TTVs. 

Another assumption considered in our study is the existence of a giant planet in a 
short-period orbit around a low-mass star. In general, the low masses of circumstellar 
disks around these stars, combined with the short lifetime of the gas in such disks, 
casts doubt in the possibility that giant planets may exist in these systems and in 
close-in orbits (Laughlin et al 2004). The fact that giant planets have been observed 
around M stars [e.g., GJ 876 with two Jupiter-like planets and a Uranus-mass body 
in approximately 30, 60, and 120 days orbits (Rivera et al. 2005, 2010), or HIP 57050 
with a Saturn-mass planet in a 42 days orbit (Haghighipour et al 2010)] suggests that 
these planets might have formed at the outer regions of the disk where more material 
was available, and migrated to their current orbits. As shown by Zhou et al (2005), Fog 
and Nelson (2005, 2006, 2007a & b, 2009), Raymond et al (2006), Mandeh et al (2007), 
and Kennedy and Kenyon (2008), terrestrial planets may form during the migration of 
a giant planet around a Sun-like star. These planets may be captured in resonance with 
the giant planet and migrate to smaller orbits. Recently Haghighipour and Rastegar 
(2011) carried out similar simulations around M stars. By allowing the giant planet 
to migrate to orbits with periods of 3-5 days, these authors showed that around low- 
mass stars, resonance capture during the migration of the giant planet seems to be 
the most viable mechanism for the formation of a system that consists of a transiting 
giant planet and a small interior perturber. In case of an outer perturber, although it 
is possible for small planets to migrate and capture larger bodies in resonance (e.g., GJ 
876 where the outermost planet is a Uranus-mass object in a 4:2:1 resonance with the 
two inner Jupiter-size planets), it is more probable that the smaller, outer perturber is 
accompanied by another larger planet in an exterior mean-motion resonance. How such 
a second exterior perturber affects the results presented in this paper is the subject of 
our currently on- going research. 
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Fig. 1 Graphs of TTVs for a three-body system with a O.367M0 star, a Jupiter-mass tran- 
siting planet in a 10-day orbit, and an Earth-mass perturber. The top and middle panels show 
the TTVs when the perturber is in a 16.11-day orbit. The bottom panel corresponds to near 
an interior 1:2 resonance. As expected, the amplitude of TTV is enhanced by an order of 
magnitude in this case. 
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Fig. 2 Graphs of TTVs for the system of figure 1 when the perturbing planet is near an exterior 
3:1 MMR. As shown here, increasing the time of integrations results in the development of 
long-term secular variations in the amplitudes of TTV. The top panel shows the results for 
250 days, the middle panel is for 1000 days, and the bottom panel is for 10^. 
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Fig. 3 The graph of the eccentricities of the transiting (gray) and perturbing (black) planets 
for the system of figure 2. 
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Fig. 4 The maximum amplitude of TTVs for the system of figure 1 when the perturber is 
near different low-order interior and exterior mean-motion resonances. 
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Fig. 5 Graph of TTVs for a Jupiter-mass transiting planet in a 10-day orbit around a O.367M0 
star in terms of the mass of the perturber. From top to bottom, Ppr/Ptr ~ 2, 2.5, 2.8, respec- 
tively. Note that for the top graph where the orbits are near a 2:1 MMR, the values of TTVs 
on the vertical axis are multiplied by 100 s. 
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Fig. 6 Graph of TTVs for a Jupiter-mass transiting planet in a 10-day orbit around a O.367M0 
star in terms of the eccentricity of the perturber. The mass of the perturber is IMq. Note that 
the ratio of periods as given in the figure are approximate values. 
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Fig. 7 Graph of the amphtude of TTV (top curve) for the system of figure 6 when the system 
is near a 2:1 MMR. 
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Fig. 8 Amplitude of TTV in terms of the eccentricity of the perturber for a system consisting 
of a O.367M0 star and a Jupiter-mass transiting planet. The perturber is Earth-mass and near 
an interior 1:2 MMR. Each line represents different time of integration. 



22 



t 1 




Ptr [days] 




Fig. 9 Amplitude of TTV for different values of the period of a Jupiter-mass transiting 
planet. The star is O.367M0. The graph on the top corresponds to Ppr/Ptr ~ 2.8, and the one 
on the bottom is for Ppr/Ptr ~ 2. The lines in each graph correspond to different values of the 
mass of the perturber (from top to bottom: 10, 6, 3, 1, 0.3 M^). 
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Fig. 10 The graph of the amphtude of TTV in terms of the mass of the transiting planet for 
two cases of near 1:2 and 2:1 MMR. The star is O.367M0 and the perturber is Earth-mass. 
The transiting planet is in a 10-day orbit. 
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Fig. 11 Top: The graph of the amphtude of TTV for the system of figure 10 and for different 
non-resonant cases, and high order resonant orbits. Note that in the case of Ppr/Ptr ~ 2.25, 
the values on the vertical axis are multiplied by 2. Bottom: Graph of the median eccentricity of 
the perturbing planet for a non-resonant system with a period-ratio of 2.8 and for the system 
in 5:2 MMR. 
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Fig. 12 Graphs of the amphtude of TTV in terms of the mass of the central star. The 
perturber is a Jupiter-mass planet in a 10-day orbit. The top panel corresponds to a near 
exterior 2:1 MMR and the bottom panel shows a non-resonant system where Ppr/Ptr ~ 2.83. 
The graphs in each panel correspond to a perturber with a mass of 10, 6, 3, 1, and 0.3 M0, 
from top to bottom. 
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Fig. 13 Graphs of the regions where a IM0 perturber creates TTVs with amphtudes of 20 s 
or larger on a Jupiter-mass planet. The mass of the central star is O.252M0 in the top panel 
and O.486M0 in the bottom. The numbers on each graph are the mean-motion resonances. 
The un- numbered curve is the transiting planet. 
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Fig. 14 Similar graph as in figure 13 for a O.367M0 central star. The top panel corresponds 
to an Earth-mass perturber, and the bottom panel is for a IOM0 super-Earth. The numbers on 
each graph are the mean-motion resonances. The un- numbered curve is the transiting planet. 
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Appendix: TTVs for a system with mstar = M0,mtr = l^j, ''^pr = l^e, Ar = lOdays. 
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Appendix: TTVs for a system with rustar = MQ,mtr = IMj, rupr = l^e, Ptr = lOdays. 
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Appendix: TTVs for a system with TTZstar = MQ^mtr = IMj, nipr = 1^©, Ptr = lOdays. 
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Appendix: TTVs for a system with mstar = M0,mtr = l^j, ''^pr = l^e, Ar = lOdays. 



